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Random advection of Lagrangian tracer scalar field 6{t, x) by a one-dimensional, spatially smooth 
and short-correlated in time velocity field is considered. Scalar fluctuations are maintained by a 
source concentrated at the integral scale L. The statistical properties of both scalar differences and 
the dissipation field are analytically determined, exploiting the dynamical formulation of the model. 

The Gaussianity known to be present at small scales for incompressible velocity fields emerges here 
at large scales {x 3> L). These scales are shown to be excited by an inverse cascade of 6^ and 
the probability distribution function (PDF) of the corresponding scalar differences to approach the 
Gaussian form, as larger and larger scales are considered. Small scales {x L) statistics is shown to 
be strongly non-Gaussian. Collapse of scaling exponents for scalar structure functions takes place: 
moments of order p > 1 scale all linearly, independently of the order p. Smooth scaling is found 
for —1 < p < 1. Tails of scalar differences PDF are exponential while, at the center, a cusped shape 
tends to develop when smaller and smaller ratios x/L are considered. The same tendency is present 
for scalar gradients PDF with respect to the inverse of the Peclet number (the pumping-to-diffusion 
scale ratio). The tails of the latter PDF are however much more extended, decaying as a stretched 
exponential of exponent 2/3, smaller than unity. This slower decay is physically associated with the 
strong fluctuations of the dynamical dissipative scale. 

PACS numbers 47.10.-tg, 47.27.-i, 05.40.-tj 


1. INTRODUCTION 

Small-scale statistics of a passive scalar advected by a large-scale incompressible velocity field is an old problem in 
turbulence theory ||^,||. Scalar fluctuations are maintained by a large-scale forcing, with typical scale L. According 
to the classical picture of direct cascade of the scalar, the injected scalar is mainly transferred downscales to the 
convective interval and then to the dissipative region. For smooth velocities, the statistics of this scalar transfer can 
be analyzed systematically and has been characterized in much detail in ||,^,|j. The core of the one-point scalar 
probability distribution functions (PDF) is Gaussian with variance O(logPe), where the Peclet number Pe is very 
large |^. Far tails of the PDF decay exponentially The physical key ingredient at the basis of these results 

is that material lines are stretched, i.e. the maximum Lyapunov exponent A for particles separation is positive. 
Typical trajectories will therefore be exponentially stretched and dynamically contracted trajectories are so rare, 
that they can only affect the extreme tails of the statistics. On the other hand, much interest has been recently 
attracted by Kraichnan model ||^ for its intermittent scaling behavior ||^,|^,^ . The picture emerging there is that 
dynamically contracted Lagrangian trajectories play a crucial role for structure functions scaling exponents C,p and, 
thus, for intermittency. The constant asymptotic behavior of C,p for large orders p found in |l^ comes, for example, 
from the most contracting possible trajectories and the value ^oo from nontrivial fluctuations of the degree of freedom 
constrained by incompressibility to still be dynamically stretched. 

Previous remarks have led us to investigate scalar transport in smooth compressible flow. The motivation is 
that compressibility might slow down Lagrangian separations and, thus, lead to nontrivial scaling and intermittency 
properties. These can be then analyzed systematically, using techniques specific for smooth velocities. Positive 
Lyapunov exponents are indeed characteristic of isotropic, solenoidal flow |^,^, but this property might be lost when 
compressible flow are considered. It is, for example, known that for compressible flow a substantial slow-down of 
large-times transport can take place (see Q). Since trapping effects are amplified when the dimensionality of space 
is low, we have focused our attention on the one-dimensional case. More specifically, we have considered the smooth 
limit of the one-dimensional (/-correlated in time model recently introduced in |^. In the absence of pumping and 
dissipation, any function of the tracer 0(t, x) (say temperature) is advected along Lagrangian trajectories and globally 
conserved in average (provided the velocity is temporarily fast or spatially smooth). Switching on the energy (d^) 
supply at the integral scale, one expects that a steady (or quasi-steady, as discussed in Section ||) distribution of the 
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scalar is attained. The main question raised here is how trapping effects due to compressibility affect the redistribution 
of energy among the scales and the intermittency properties of the scalar field at the stationary state. To answer this 
question, we exploit the dynamical formulation of the model to calculate the statistical properties of 9 both up and 
downscales, i.e. at the scales smaller and larger than the integral scale L. Since the scalar is a tracer in the velocity 
field and the velocity is smooth, the problem reduces to studying Lagrangian separations statistics. 

The major physical difference appearing with respect to the incompressible case is that the maximum Lyapunov 
exponent for Lagrangian separations is negative. This means that, along typical trajectories, distances are mostly 
contracted and the stretching process is strongly intermittent in time. As in [^, the second and higher powers of 
the distance R(t) between Lagrangian trajectories grow exponentially, while its low-order positive moments decay 
exponentially. This is the dynamical origin of the major results found in this paper: inverse cascade and Gaussianity 
at large scales and extreme intermittency at small scales. 

Scalar correlations are indeed essentially governed by the statistics of the time spent by particles at distances smaller 
than L. Consider then two particles initially separated by a distance x ^ L. Despite of the fact that x ^ L initially, 
the distance R{t) will typically reduce to 0{L) by a large time ~ log (x/L) /|A|. The consequence is that even scales 
much larger than L are strongly excited and this is the dynamical hint of the inverse energy cascade. Moreover, for 
moments of order n <gC log (x/L), relative fluctuations around the previous typical time are small and this leads to 
the Gaussianity, e.g. of the core of scalar differences PDF. On the contrary, the statistics at small scales x <C L is 
associated with the stretching process. The time to reach the integral scale L strongly fluctuates and this is reflected 
in the intermittency of both scalar differences and gradients. 

Note that the inverse cascade of the scalar taking place here differs in one important respect from other known 
examples of inverse cascades (say, energy cascade in 2d Navier-Stokes turbulence @ and number of waves cascade in 
wave turbulence |p^ ) : no flux of another integral of motion (like enstrophy or wave energy) is present. The origin of 
the inverse cascade found here is purely dynamical and associated with trapping effects. An interesting consequence 
of the inverse cascade is that the equation for velocity differences PDF, derived here exactly from the dynamics, 
coincides with the one without dissipative anomaly (operator product expansion, which may result in the anomaly, 
was proposed in [Q in the context of the Burgers turbulence; see also |19 for possible extensions to the passive scalar 
turbulence). The absence of anomaly is indeed a consequence of the inverse cascade, preventing the rare trajectories 
emerging from the dissipative range to affect the convective range behavior. 

Strong downscales intermittency emerges all over the quantities calculated in Sections IV-VI. Moments of scalar 
differences of order n > 1 all scale linearly with x, independently of the order. This collapse of exponents carries over 
to the dissipation field e, that has all its integer moments scaling with the same power of the Peclet number. Smooth 
scaling is observed for low-order moments of both scalar differences and dissipation. Very large fluctuations of these two 
quantities behave, however, quite differently. Scalar differences PDF has indeed a Lorentzian shape for values smaller 
than unity and exponential tails. The tails of the dissipation field PDF are, on the contrary, stretched exponentials 
with exponent 1/3 (and not 1/2). The additional probability for these strong events comes from fluctuations of the 
dynamical dissipative scale. Gomparison of the n-th moment of the dissipation field with the 2n-th moment of scalar 
differences establishes indeed the effective viscous scale. This appears to be a strongly fluctuating quantity, growing 
factorially with n EJ. This factorial dependence is the cause of the 1/3 stretched exponential. 

Even though the explicit calculations are quite lengthy, the underlying technical ideas, that make the analytical 
calculations doable, are simple to explain. Since the scalar is passive, its n-th order correlation function can be 
presented in terms of a matrix element of an auxiliary “Quantum” mechanics. There are then two import steps in the 
evaluation of such matrix elements (and thus of the scalar PDF’s). First, smoothness of the velocity field (stretchings 
and compressions are uniform on all the particles) allows to reduce (Section ^ and Appendix A) the multi-particle 
to a one-parameter problem. The fact that the Lyapunov exponent is negative clearly emerges in this procedure. 
This leads, after the very direct calculations of Sections ^ and |^, to a compact expression for the PDF of the scalar 
differences in the convective interval. Here, the temporal dynamics of the fluctuating degree of freedom is local, while 
the locality is lost in the dissipative range. The second step comes then into play. To describe the convection-diffusion 
interplay in Section and Appendix B, we use a scale separation procedure. The temporarily nonlocal (diffusive) 
and local (convective) dynamics of the fluctuating degree of freedom are well separated by time to, 1 to ^ ln[Pe], 
if the Peclet number Pe is large. Independence of the resulted average (say, the gradient’s PDF) over both local and 
nonlocal domains on to, and smallness of the neglected terms with respect to the inverse Peclet number, justify the 


^An essential enhancement of the dissipative scale was observed also in where the tail of the scalar pdf was studied for 
the incompressible case by an instanton technique. 
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scale separation procedure. 

The plan of the paper is as follows. In the following section, the one-dimensional passive scalar model is briefly 
recalled, its relevant time and length-scales are discussed, and the inverse cascade issue is explained from consideration 
of scalar pair correlation function. The dynamical formulation associated with the passive scalar equation is the subject 
of Section and the Appendix A. The latter is based on the Martin-Siggia-Rose formalism, while the former is in 
terms of particles formalism. A key point arising in both procedures is that compressibility couples the dynamics with 
a global mode. This mode must then be taken into account in order to get the dynamical formulations. Multi-point 
correlation functions of scalar gradients are discussed in Section |^. Scalar differences PDF is described in Section ||, 
where behaviors up and downscale with respect to L are considered in two different Subsections. To describe the 
advection-diffusion interplay, we develop a scale separation procedure in Section which is devoted to the PDF of 
dissipation. We use the scale-separation procedure also in Appendix B to study the question, how does the steady 
regime for the pair correlation function of the gradients (discussed in Section ^ form ? The final section is reserved 
for conclusions and discussion of questions which may be of a general relevance for other problems in turbulence and 
physics of disordered systems. 


2. THE MODEL. PAIR CORRELATION FUNCTION AND INVERSE CASCADE. 

Our aim here will be, first, to briefly recall the equations of the one-dimensional model introduced in Q, and then, 
solving the equation for scalar pair correlation function, show the effects of compressibility on the redistribution of 
energy among the scales and how does the inverse cascade of the scalar work. j-. 

The advection-diffusion equation governing the evolution of the passive scalar 9{t, x) is 0 

dtO + ud^e = Kdle + f. ( 1 ) 

The velocity field u and the force / are both assumed to be Gaussian and ^-correlated in time. The force has 
correlation function 


x) fit', x')) = 6it- t')x ^ , (2) 

spatially concentrated at the integral scale L. The velocity has zero average and correlation function 

{u{t,x)u{t',x')) = \Vo — Six — x')]5{t — t') with Six) = . (3) 

The specific smooth case considered here corresponds to 7 = 0. The scaling behavior of the structure function S 
persists up to the infrared cut-off Lu, the largest scale in our problem. The scale independent part of the velocity 
correlation function Vq is estimated by the infrared cutoff of the velocity field squared, ^ L^. For scales larger than 
Lu, velocity correlations decay to zero, i.e. the structure function saturates to the constant value Vq- The equation 
of motion for the gradient field ui = d^O immediately follows from (|^) 

dtU} + da: iu uj) = ndlio + dxf. (4) 

The d-correlation in time of both the velocity and the forcing allows to derive closed equations of motion for equal¬ 
time correlation functions [||. It is for example easy, e.g. using Gaussian integration by parts (see [ 0 ) , to derive the 
equation for the second-order correlations Llit,x) = (a;(t, x)a;(t, 0)) and Fit,x) = (d(t, x)0(<, 0)) 

dtVt = dl[iSix) + 2K)n\-dlx, (5) 

dtF = iSix)+2K)d'^F+ x- ( 6 ) 


^Note that there are two types of passive fields for compressible flow. Lagrangian tracers, like entropy or temperature 
(provided pressure is slowly varying in space and time), are conserved along Lagrangian trajectories in the absence of diffusion 
and pumping. The local maxima of the field do not grow in the absence of pumping. Concentration fields, e.g. of a pollutant, 
are, on the contrary, only globally conserved and their maxima can be amplified. The equations for the two types of fields 
differ by the position of the space-derivative in the velocity term. In our one-dimensional case, the two possibilities correspond 
to the 6 and the u> fields, respectively. 
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It is convenient to consider first equation (|^) and then recover the correlations of the scalar by integration. From the 
very definition of the correlation function, it follows that the solution of i) is even in x. Looking for a stationary 
solution, there is another boundary condition needed to fix the integration constants. This comes from the dynamics. 
Let us indeed consider the situation when the system starts from rest (0 = 0). At any subsequent time, the solution 
satisfies f fi dx = 0, where the integral is taken over all the space. Note that such a condition is consistent and actually 
dictated by the dynamics: the integral of the correlation function of the gradients is conserved in the evolution on 
account of the double space derivative in the r.h.s. of (||). It is then easy to find the stationary solution of 


n(a:) 


X(x) _ C 

2k + S(x) 2k + S(x) 


where 


Jx/i^K + S) 
f 1/ (2k + S) 


(7) 


(see Appendix B for a dynamical derivation of (|^)). The expression (^ is very illuminating for several reasons. Let 
us first consider the case where no infrared cutoff is present. The integral f \/(2k + S) is then convergent and, 
for small molecular diffusivity, varies as Xj^[ k. For the sake of concreteness, let us specifically consider here a forcing 
that is regular at the origin. The dominant terms for distances much smaller than the integral scale L are then 


fl{x) 


lx"(o)l 

2{2k + S) 


[aL^/K — x"^] , 


( 8 ) 


where a is constant 0(1), dependent on the detailed form of the pumping. The first and the second term dominate 
respectively at scales smaller and larger than L/v^Pe, where the Peclet number Pe = is supposed to be very 

large. Note that this scale is still much larger than the dissipative scale L/Pe. The most interesting aspect of (^) 
is that the dissipation Kr2(0) vanishes as 1/Pe, i.e. there is no direct cascade. The energyH is actually transferred 
upscales by an inverse cascade, as also emerges from the behavior of scalar correlations. Let us indeed insert into (||) 
the correlation {9(t,xi)0(t^X2)) = S-^oo ~ V^) with the w-correlation function having the expression (0). It 
is easily checked that the energy satisfies dt{9‘^) = O, i.e. it grows linearly with time. Note that the growing mode 
is constant in space and thus disappears when differences or gradients are considered. The effect of the advective 
term udx9 is therefore to transfer energy upwards in the scales. Since dissipation is quadratic in the wavenumber, the 
energy on the large scales cannot practically be dissipated and it is piled-up. This is at the origin of the linear growth 
in time of ( 0 ^). C 7 ^ 0 correponds to the inverse cascade, which therefore holds generically for any scaling exponent 
7 < 1 in (H) (if the source function y is not exceptional). 

Let us now introduce an infrared cutoff Lu in the velocity field. Since S saturates to a constant, it is clear now 
that the integral f 1/(2k + S) diverges. The constant C in must then vanish. In the presence of a finite cutoff 
Lu) there will then be a very long intermediate in time asymptotic where for scales the behavior without cutoff 

is observed. However, after a very large time ^ (L„//c)^ the dynamics changes: the inverse cascade stops, (0^) 
saturates to a finite value, the system starts dissipating a finite amount of energy in the limit of large Pe and the 
correlation function U,(x) tends to the solution ( 0 ) with (7 = 0. The finite contribution — J x/(2«: + S) needed to 
ensure the zero integral condition comes from a strip of negative values that becomes more and more extended with 
time and whose amplitude tends to vanish. More details on the dynamics of the pair correlation function of the scalar 
gradients at infinite times in the presence of an infrared cutoff Lu may be found in Appendix B. 


3. DYNAMICS IN PARTICLE (LAGRANGIAN) FORMALISM 

In this Section we shall discuss the dynamical formulation of the equations of motion. The goal is the same as 
in reduce the calculation of simultaneous scalar statistics to averaging of functionals of the random-in-time 

strain-vorticity matrix. This reduction is crucially based on the fact that the velocity field is smooth (7 = 0) and can 
be performed for any smooth flow, independently of its compressibility and the dimension of space. In the specific Id 
case, no matrices are obviously involved and one is left with averaging of a single scalar field. Compressibility makes 
however the derivation slightly more involved and some care must be taken in the ordering of the advective term. This 
emerges, in particular, in the non-vanishing average of a in the weight (|^) for the Lagrangian trajectories (^). The 


^It is important to mention that, in the absence of pumping and dissipation, the average over the velocity of any function of 
the scalar {f{9)) is conserved, although the average over all the space of f{0) itself is not conserved in the particular realizations. 
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dynamical formulation is derived here using particle formalism. The equivalent derivation using Martin-Siggia-Rose 
field formalism is reserved for the Appendix. 

Equation (|^) for the passive scalar can be presented in the form 


e[T-x) 


Texp 


V{t'] x)dt' 


1 

x)dt = J dt J dy^{t,T\x,y)(j){t\y), 


(9) 


where it was supposed that no pumping was supplied at negative times. In (^), the operator 7^(t; x) = —u(t; x)da:+Kdl, 
time-ordered exponential is denoted by T exp and the function dt can be expressed by using Lagrangian trajectories 
as 


p(T)=x 

'^{t,T;y,x)= J Vpexp 

p(t)=y 



/4k 


( 10 ) 


This formula expresses the simple fact that Lagrangian trajectories are fixed by the velocity u and smoothed by the 
molecular diffusivity k. One can express in the more convenient Hamiltonian form 


^N-l N-l 


'i/{t,T;x,y) = / tt tt dpnexp 


n—0 n—1 


N r 


AE 

n—1 


A 


Vn—l^Pn Pn—l) Pn—l'^{in: Pn^ 4“ 


Pn= p{tn=t + nA), p„=p{t + nA), A = 


T-t 

N-l 


( 11 ) 

( 12 ) 


where the momenta integrations (dpk) run along the imaginary axis and regularizations have been specified. Using 
the property that both velocity and pumping are Gaussian and have correlation functions (^) and (^, we can easily 
perform the averages in the 2n-th simultaneous product of the scalar field 9{T; x). We thus obtain 


F{T;Xi,--- ,X2n) = {0{T,Xi) ■ ■ ■6{T,X2n)) 


/ p 2n 

di / It dy{R.{T - t; Xi, yi) [F{t\yi, - , y 2 n- 2 )x{y 2 n - y 2 n-i) + permutations], 


(13) 


where the eddy-diffusivity resolvent is defined as 


2n j \ 

TZiT;x^,yi) = Jl/ ^{T,0;Xi,yi) 
i=l \ I 

Pi(T)=yi 


Pi(0)=Xi 


'Dpt{t)'Dp^{t) exp I J dt 


Pik - -jPi{Pi - PjfPj + Kpf 


(14) 


Here, one inverses the direction of time in comparison with & { t ^ T — t) and thus the convective term (p^p^) is 
regularized in a way such that its p-dependent part is retarded with respect to the p-dependent one. The Hubbard- 
Stratonovich transformation of the diffusive term gives 


Pi{T)=Vi 


'R-{T]Xi,yi) = J Vp^{t)Vpt{t)'D^i{t)exp[-S - S^] 


T 

S = J dt 
0 


Pi{0)=Xi 


Pi pi Pj Pi^i 




4k 


e,dt. 


(15) 


(16) 


In the integration in (13) of the resolvent, the correlation function F and the pumping correlation y appear. Both 
depend on the difference of the coordinates only. This means that one can simply integrate with respect to collective 
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variables, say p = ^Pi and the collective momentum P = Pi/{‘2n). To get rid of the superfluous degrees of free¬ 
dom, let us come from the old set of variables {pi, • • •, p 2 n',Pi, ■ ■ ■ ,P 2 n} to the new one {p, p 2 , - ■ ■, P 2 n', P,P 2 , ■ ■ ■ ,P 2 n}- 
The new momenta Pi = Pi — P are considered in the system commoving with the ” center of mass” and the positions 
Pi = Pi — Pi are with respect to one, e.g. pi, taken as reference. The action S can then be decomposed as 5 = Scoi + S, 
where 



-Pp-P^e. 

2=1 


(17) 

(18) 


2n 

Let us now consider the integral f fl dyiTZ{T ; Xi, yi)f{yi) where, as in ([l^, / is a function of differences only f{yi) = 

i=l 

fiVi — y)- Collective degrees of freedom p and P are easily integrated. The principal point for this integration is the 

2n 

absence of a “potential” p dependence in the action S. The integral is then reduced to / fl dyiR{T ; ii, yi)f{yi + Pi), 

i>l 

where the effective resolvent 


Pi{T)=yi 

'R{T]Xi,y{)= J 'Dp{t)'D^i{t)'Da{t)exp[-Sa - S^]Y\_S 

Pi(0)=5i 

1 2 


(m) _ (m— 1 ) \ 

Pi - p -( 19 ) 


and the measure of averaging = '22m /4. In order to derive (^^, we have decomposed the quadratic overpp 

term by means of Hubbard-Stratonovich trick, introducing an additional collective integration over d. The integration 
with respect to momentum pi is already performed in the last line of (H) (the effective paction appears to be linear 
in p in the result of the Hubbard-Stratonovich transformation). The continuous versiongj of the equation under the 
i5-function sign in dil) and the measure of averaging over a are 


pi — ^Pi Cl: ^ — d 1, 

T 

5 ,^ = i J dt [a + if . 


( 20 ) 

( 21 ) 


To see the relation between ( |19[) and (^|^) one can check, particularly, that both discrete and continuous versions 
give {pi{t)) = Pi(0). The formal solution of (^) is 


( 22 ) 


It is finally easy to recalculate the 2n-particle correlation function of the scalar from the 2n-particle resolvent 


t 

■ t 

[ dPW-\t'){C^-^l), w{t) = exp 


J 

0 

J 

.0 


F{T-xi,- ■ ■ ,X 2 n) = ^^[T; {cr(t)};a:i - X 2 ] ■ ■ ■ :i[T; {a{t)};x 2 n-i - X 2 n] + permutationsy 


l[T;{a{t)};Xi- Xj]= / dt( x 


1 

!■ 


w(t p +w{t) [ 


(23) 

(24) 


_ 


where averages over a{t) and ^i(t) are fixed by the measures exp[—5^] and exp[—5^] defined in (^ and (0), 
respectively. 


'^By continuous version we understand, particularly, a symmetrical smearing of the temporal 5-function in the cr-field correlation 
for a small but finite (which is still much larger than the temporal slice A) width. 
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4. CORRELATION FUNCTIONS OF SCALAR GRADIENTS IN THE CONVECTIVE INTERVAL 


Using the results of the previous Section, the dynamical expression of correlation functions of the scalar gradient 
u){t,x) = dxd{t,x) can be simply found differentiating with respect to all spatial arguments. Here, we shall be 
interested in the behavior of these simultaneous correlation functions for distances such that molecular diffusivity 
can be neglected. We first derive a general formula valid for arbitrary form of forcing correlation and then treat 
more specifically the case of an exponential pumping. The resulting expression shows that the ratio between the 
irreducible and the reducible contributions to the 2n-th correlation function grows as (L/a;)" This evidences 
both the non-Gaussianity of the field and the fact that it increases going downscales in the convective interval. The 
same phenomenon takes place for generic forms of the pumping. The discussion of the range of scales where these 
considerations hold, is postponed to the end of the Section. 

It follows from that scalar gradients correlation functions are given by 


{u;(T,Xi) ■ ■ • W(T, X2n)} 

(-ir 

L2n 


■ T 

I dtW\t)x" 



■ T 

[ dtW^t)^' 

J 

.0 

L J 


J L 

.0 


W{t) 


^2n ^2n—1 


permutations, 


(25) 


where molecular diffusion effects have been neglected. The only averaging left in (25) is then with respect to the a 
statistics defined by (^^. A possible way of performing this average is to introduce the auxiliary object 


A(sn2, • • • , S2n,2n-l) = ( exp 


W^it) Si,2X' 


W{t) 


Xi — X2 


+ • • • + S2n,2n-lX 


W{t) 


^2n ^2n—l 


dt 


). (26) 


Differentiating (|26D over s variables, (|^) is clearly reproduced up to the sign and the L dependent prefactor. Inserting 
the weight (|^) into (| 2 ^), one can easily recognize the path-integral structure associated with the time evolution of 
the “quantum mechanical” Schrodinger equation, having Hamiltonian H = —d^ — C/(exp 77 ) exp(2/;). The “potential” 
U appearing in the Hamiltonian is 


Uiy) = si.2x" 



+ • • • + S2n,2n-lX'' 


^2n 

y — 


^ 2 n —1 


(27) 


t 

and the “space” variable 77 is defined as 77 = J dt'a{t'). Using standard notation for quantum mechanics matrix 

0 

elements, it is easy to check that (@ can be presented as 


A = exp[—T/ 4 ]^i 5 (r 7 ) 


-TH 


,-v/2 \ _ 


e-^/4n(r;77) 


r}—0 


(28) 


where the ’’wave-function” n(r; 77) satisfies (dr — i7)n = 0, and the initial condition is n(0; 77) = We can now 

remark that the potential part of the Hamiltonian vanishes at 77 —*■ —00, while the initial condition does not. The 
resulting asymptotic behavior at large times T will then be 


n(T;r 7 = ln[ 7 /]) 


„T/4 


Hn 




[dl + U{y)] Ho (7/) = 0. 


(29) 


The new variable y = exp 77 has been introduced. The boundary conditions for the spatial part Hq are easily derived 
from those for H: it should tend to unity for y 0 and the ratio Hg /^yy should vanish for y tending to infinity. It 
follows from (p^) that the auxiliary object A is simply the function Hg calculated at 77 = 1. 

Derivatives of A at s = 0 are needed for the calculation of gradients correlations (p5|). It is then convenient to 
present the solution of (p^ in the form of an expansion with respect to the “potential” U: 


A(si,2, • • • , S2n,2n-l 



dy2kU[y2k]- 


(30) 
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Only the n-th order term of the expansion in (30) actually contributes to the 2n-th order scalar gradients correlation 
function. Its final expression reads 

1 OO 


(— 1 )” , f f 

(Wi • • • W2n) = ^ j^2n X! / / ^^2%" 


y2- 


^ki ^k2 


y 2 n -2 OO 

X j dy 2 n-l J dy 2 nX 

0 y2n-l 


0 yi 

^k2n ~ ^k2n-l 


V 2 OO 

J dya J dyix" 

0 ya 


Vi 




y2n 


(31) 


where summation is performed over all the splittings of the set {fci, • • •, fc 2 „} into n ordered pairs. 

Formula ( |M| ) holds for any pumping x- Let us now specifically consider the case of the exponential pumping 
Xe{z) = exp[—|z|], where the integrals appearing in m can be easily performed. For the second-order correlation, 
we obtain that the dominant contribution for small distances X 12 /L is — l/(a;i 2 L), in agreement with the solution ( 0 )- 
For the fourth-order correlation function, we find 


{uj{xi,t)uj{X2,t)uj(x3,t)uj{x4,t))e = Fe(xi2,X34) + Fe(xi3,X24) + Fe{Xi4, X 23 ), Ff,{x,y) ~ 


(32) 


where Xij = 


xy{x + y)L' 

Xi — Xj\, the subscript is intended to remind that this specific expression holds for the exponential 
pumping and only dominant terms in XijjL have been retained. Distinguishing between reducible (Gaussian) and 
irreducible contributions into (32), one observes that the irreducible part is L/a; ^ 1 times larger. More generally, 
if X denotes the typical distance among the various particles, i.e. Xij ^ x, then (cci • • - 1 x 1211 ) 1/ {x‘^^~^L). Both 

the correlation functions and the degree of non-Gaussianity of the scalar gradient (ratio of the 2n-th moment to its 
reducible part) are then growing with Xij going downscales. 

Let us finally come to the range of validity of the previous convective arguments. The exponential form of the 
pumping is a special one, being not regular at the origin. The first term of its expansion at small distances is linear 
and not quadratic. This affects the dependence on x of the second-order correlation. For a regular pumping, the 
dominant term would indeed be constant, as follows from (||). On the contrary, one can check using (|^) that the 
dependence of correlations of order > 4 on a; and L remains the same as for the exponential pumping. The regularity 
at the origin of the pumping also enters the range of scales where neglecting molecular diffusivity effects is allowed. 
Performing the small distance expansion, as in for a regular pumping, a logarithmic correction oc logPe appears. 
The ultraviolet criterium of applicability of the previous convective considerations for the exponential pumping is 
then X ^ -v/Klnpe]. For a pumping regular at the origin, the criterium is the same as in (^, i.e. x ^ 


5. SCALING AND PDF OF SCALAR DIFFERENCES 


Scalar structure functions S 2 n{x) = {{0{T,x) — 0(T, 0))2") can be easily expressed in terms of scalar correlation 
functions F by taking the appropriate combinations of them. A dynamical expression for S 2 n can thus be derived 
directly from (^^, obtained in Section ^ for the F’s. This is however not a very practical procedure. Each of the 
contributions, appearing in the sum giving the structure functions, contains indeed the constant mode. As it was 
discussed in Section ^ and as it also appears from the dynamical expression ( ^ ) , this mode grows linearly with the 
observation time T. It is just in the whole sum that these divergent contributions are canceled, thus leaving the time- 
independent final result for structure functions. It is then more convenient to restore structure functions directly from 

X X 

scalar gradients correlations as S 2 n{x) = f dxi ■ • ■ f dx 2 n(‘ 2 Ji ■ • ■ 0 J 2 n}- An important question is whether, to calculate 

0 0 

S 2 n in the convective interval, we can avoid to take dissipation explicitly into account or not. This means essentially 
asking whether it is enough to know just the convective expressions for gradients correlations or their whole behavior 
is needed. This point can be tested by simply taking the convective expressions for (oji • ■ ■ i 02 n) found in the previous 
Section and inserting them into the integral expression for S 2 n- One can then check that all the integrals for any 
structure function are convergent on the ultraviolet and dominated by the infrared side of the convective range. 

The expression for structure functions in the convective interval is then 


S 2 n{x) = 

(2 n-l)!!n!2" J dyi J dy 2 


xj L 00 2/2 00 y2n-2 OO 

J ay, I ... i ay,„., J . (33, 

0 yi 0 ya 0 yan-i 


Vln 


















where we have already performed the 2n integrals over the dxi’s. The whole set of equations ( |3^ can be recast into 
the more compact equation 



x(0)-x(f)])2(|,A) = 0, 


(34) 


for the generating function Z(x,X) = {exp {—iX69x)) of scalar differences S6x taken at the scale x. From this very 
definition, it follows that the function Z must tend to unity for vanishing x and, for the convergence of the integrals 
in (33), it should grow slower than linearly at infinity. 

It is worth to remind that (|^) was found as the result of an accurate dynamical procedure: we first calculated 
correlation functions of the scalar gradient for all points being separated; the resulting gradients correlation functions 
were then integrated to obtain structure functions and, finally, the generating function for scalar differences was 
reconstructed from its moments. We thus avoided to explicitly handle diffusion, paying for this the price of taking 
many particles into consideration. The closed differential equation (|^) for scalar differences generating function 
emerges as the result of this procedure. On the other hand, one could generally (also for the non-smooth case, i.e. 
7 ^ 0 in (H)) derive the following unclosed Fokker-Planck equation for the generating function 


{x^-^dl - X^ [x(0) - xix/L)]) Z{x/L, X) = J [dlOi - dle^] exp [A [01 - 02 ]] 


(35) 


This equation is simply obtained averaging the equation of motion (^ for the scalar at two reference points. The 
smooth (7 = 0) limit of (^) differs from ( ^^ by the r.h.s. dissipative term. There is a general expectation that 
this term may remain finite even in the limit k —> 0 , thus providing a nonvanishing anomaly in the terminology of 
Polyakov |^. Equation (^), derived microscopically without any conjecture, shows then the absence of anomaly for 
the one dimensional smooth flowO. Physically, the absence of anomaly is associated with the vanishing direct flux of 
0^ in the limit of infinite Peclet numbers. This point, already appeared in Section |^, will emerge even more clearly 
in the full analysis of the dissipation field in the next Section. In the two following Subsections, solutions of ( ^^ at 
X ^ L and a; <C T will be discussed. 


A. Upscale interval 


Let us first consider scales larger than the integral scale L. The asymptotic solution for the generating function, 
which can also be found from (^ replacing x( 0 ) — xi^jL) by x( 0 ), is 


/xy/2-Vl+4x[0]AV2 

Z=(-) . (36) 

The inverse Fourier transform of (|^), calculated in the saddle-point manner (the large parameter is x/L or, equiva¬ 
lently, large values of l(50a;l), gives 




1 

2 v'^x[ 0 ] ln[a;/L] 


exp [-j/2/(4x[0]ln[a;/L])] , 


exp 


-|y|/(2v^) 


jyj < \n[x/L], 
\y\ > ln[a;/L]. 


(37) 


Thus, the arguments for Gaussianity presented in for the incompressible case can indeed be reversed and applied 
here for the upper interval. Structure functions of orders much less than ln[a:/L] scale indeed logarithmically, the core 
of scalar differences PDF is Gaussian and the PDF’s tail is exponential. 


B. Downscale interval 

We shall now obtain a general formula for the PDF at ?/ <C 1, no matter how y and xjL relate to each other, 
provided both of them are small. Indeed, replacing the x dependent “potential” by the first term of the expansion 


®We acknowledge A. Polyakov for attracting out attention to the matter discussed iu the paragraph. 
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over xjL, and solving the resulting differential equation with the same boundary conditions as before, one gets 
the simple exponential form for the generating function ZixjL^X) = exp[—A-\/x"[0]/2 x/L]. The inverse Fourier 
transform produces the following Lorentzian expression for the PDF 


p(«.)(„) = if 


TT X y'^ + |x"[0]| x^ 


(38) 


which is thus generally valid at \y\ <C 1. It results from (|38| ) that the PDF is smooth in a small region around the origin 
?/ = 0, where it can actually be expanded in y'^Lp/x^. This region extends approximately up to xjL, where the second 
behavior in 1/j/^ sets in. Note that the concavity of the PDF in this second region is upwards and remain upwards 
up to the small values xjL. It is just for very small values \y\ ^ xjL that the concavity is reversed downwards. An 
“experimental” histogram of such a PDF would then look strongly cusped at the origin. 

The tail of the PDF matching the Lorentzian (^) at |j/| ~ 1 is exponential. To see this, and generally to obtain an 
explicit analytic solution for the PDF in the whole domain of x and 59x, let us consider the following specific form of 
pumping 


l + {x/Lf 


(39) 


Making in (^) the change of variable cot[</j] = x/L, the solution of the resulting equation can be expressed in terms 
of associated Legendre functions as 


Z*(cot[(/j]; A) 


2 ^-i/ 2 (^ + 1)F^ [{u + l)/2] (cosM) 

V^TT sin[(p] 


(40) 


where the upper index i/ = —1/2 + 1/2 y/1 + 4A^. The choice of the sign for the square root is such as to ensure that 
the generating function Z{x,\) grows at infinity slower than linearly. The notation Z* in ( ^0| ) is intended to stress 
that this explicit solution was obtained for the pumping (^). Note, that (^) is particularly applicable for the upper 
interval discussed in the previous Subsection. Considering (4^ at <C 1 and A 3> 1, one indeed recovers (|3^). Using 
the integral representation (8.714) in and the doubling formula for the gamma function, ( p0| ) can be presented in 
the integral form 


Z*(cot[(/j]; A) = 


2r[(i^ + 3)/2] 

/^sin[(/5]r[(^ + 2)/2] 


cos[t] 


COs[t] — COs[(/7] 
sin[(^] 


dt. 


(41) 


The integral representation (41) is useful, since it clearly shows the analytic structure of Z* with respect to A. One 


observes, particularly, that Z* is analytic in the whole upper semi-plane, except on the line of imaginary A from 
i/2 to ioo. The path on the real axis in the inverse Fourier transform can be then deformed into the following one 
surrounding the cut: the left branch goes from ioo — 0+to i/2 — 0+, bypass of A = i/2 from below and the right 
branch goes from i/2 -|- 0+ to ioo -I- 0+. The final result is 


><“•>(!/) = T/ 


qdq 




exp 




g2) 


^(cotM;g), 


(42) 


where the function —iQ is the difference between Z* on the right and the left of the cut, i.e. v —1/2 ± ig/2, 
respectively. The general expression for Q can be derived from the integral representation (^) as 


Q{coi[ip\\q) = 


4 J cos[t]dt 

rp(WH) 

sin[(/5] 

-iq/l' 

Y'TTsin);^] J Y'cos[t] — cos[(/5] 

p(3^) 

_cos[t]—C0s[(/3]_ 

. 


(43) 


The calculation of the PDF of scalar differences in the convective interval is now reduced to the evaluation of the 
asymptotic behavior of Q in ( ^^ and then to perform the integral (^). In the convective interval, xjL = cot[</j] 1, 
i.e. the angle Lp is very close to 7r/2. The dominant expression of Q in this region can be obtained by simply expanding 
directly in (E^. Inserting this expansion into (42), we obtain 
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’ qyjl + sinh[7rg/2] exp + 


dq 


I r [(3 + zg) /4] 14 cosh [nq /2] 


Iss. 1 

TT L ^ ’ 

:exp[-|y|/2], 


1 > |y| » xjL, 

\y\ » 1 , 


( 44 ) 


where, at y ^ 1, the prefactor algebraic in y has not been considered. Varying the pumping function x(a;/L), one 
can change the number behind |i/| in the exponential, but the exponential behavior itself will never change. 

One may check that the PDF’s asymptotic for the smallest values \y\ xjL derived from is consistent with 

the general formula (^). Indeed, the respective large q asymptotic of Q is 


g{cot[(p]-q) - 2 


cos[t] 


7rsin[(/5] J y'cos[t] — cos[' 


sin ^ In 


sin ( 


cos[t] — cos[y)] 


(45) 


Substituting ( ^ ) into ( ^ ) and performing first the integral in q and then the one in t, we obtain at \y\ <C xjL 1 

1 1 cosl(|d( 1 1 i 1 _ 11 ^ V -^ I ^ I 1 _ 1 (46) 




^/2tt J y^cos[t] — cos[' 


cos[t] — COs[(/3] 


8 In^ 


cos[t] —cos[c/3] 


1 L 

TT X 


Expressions (Q) and ( ^^ are clearly in agreement with (|^), valid for an arbitrary form of pumping. 

Let us now derive the scaling behavior of scalar structure functions (|0i — ^ 2 !“) at x <C L. The scaling for orders 
a > 1 is dominated by the behavior of the PDF at values of order unity. On the contrary, for a < 1, the region in 
dominates. The resulting scaling behavior of structure functions is 


( 101 - 021 “) 


f, for a > 1; 

(f)“, for-l<o<l. 


(47) 


For a < — 1 the moments do not exist at all, since the PDF is finite at the origin. Remark that (HI) is exactly the 
same scaling as for velocity structure functions in Burgers equation. In the present model the collapse of high-order 
exponents is associated with the uniformity in space of stretchings and compressions. 


6. PDF OF SCALAR DISSIPATION AND GRADIENTS 


The n-th order moment of the dissipation field e = k {dxO{t; x))^ can be obtained from the dynamical expression 


= (2n- 1 )!!k” ( [dx^dxj^lT; Xi - Xj]]l^^^ ^ ) = (2n- 1)!! ( Q[T;{a{t)}Y 


The average with respect to the noises is easily performed in Fourier space, thus obtaining 



(48) 


(49) 


where W(t) has been defined in ( p^ and Xq is the Fourier transform of the pumping correlation function x(a;). Let 
us remind that the Peclet number is defined as the pumping-to-diffusive scale ratio Pe = L/-\/2 k and is supposed to 
be large. To obtain the stationary value (e"), the limit T ^ 00 should be considered. 

The average over a with the Gaussian weight (p^ has to be calculated in ( p^ . This is very hard to do explicitly in 
the general case. We can however exploit the presence of the large parameter Pe to develop an asymptotic theory that 
captures the dominant terms in ( p^ with respect to Pe. The important point is that, when Pe is large, there are two 
very different time scales in the dynamics. For the Lagrangian trajectories relevant to (p^, particles start very close. 
The additive molecular noise term in the Langevin equation for Lagrangian trajectories is dominant at these distances 
and remains dominant as long as the particles do not separate by a distance comparable to the dissipative scale. This 
phase of the dynamics corresponds to the formation of the integral in the square brackets in (^^ and takes place 
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on times of order unity (not scaling with Pe). Once the particles have reached the dissipative scale and enter into 
the convective region, random multiplicative effects due to the velocity become dominant. Due to the multiplicative 
nature of the dynamics, the time to go from the dissipative scale to the integral scale varies as InPe. This phase is 
associated with the growth of the terms in (|4^. For large Peclet numbers, the two processes, formation of the 
integral in the square brackets and growth of bF^erms in (^), are well separated in time. Let us then consider a 
time to satisfying 1 ^ tg ^ ln[Pe]. At the dominant order in Peclet, Q can be approximated as 


Q[T-:W{t)}] 


T 

r ^ 1 


exp 

2 f a^{t')dt' 

^ fdte^v 

9 / 

G 


to 

Pe^J P 

1 JtlL 


Pe^ 

^0 

- *0 





(50) 


where a and /3 are defined as 



*0 

*0 

p 

t' 

p 


^0 

p 

a = exp 

2 / a{s)ds 

/ dt'exp 

-2 / a{t")dt" 

, (3 = exp 

2 / a(t')dt' 


J 

. 0 

J 

0 

J 

0 


J 

. 0 


In order to obtain ( |^ ) from the original expression (^), we have made the following two steps : (5[to; {cr(t)}] has been 
neglected and the upper bound t in the integral over dt' for a has been replaced by to- Both steps are motivated by 
the time scales separation at large Peclet numbers. More precise conditions of validity of the approximation will be 
discussed later in the Section. The moments (Q") can now be obtained by taking the n-th power of (^) and averaging. 
The average over a(t) is decomposed in a small times part I?cr< = n da(t) and a large times part I?ct> = ])[ da(t). 

t<to t>to 

The corresponding weights are simply obtained decomposing the integral over t in (^ as 5< = 1/4/ [tT< + 1]^ and 

0 

in 5> the integration runs from tg to T. The great advantage of (^) is that, in the large times averaging, both a and 
/3 are just external parameters, depending neither on time t nor on cr>. Once the average over (T> is performed, we 
are then left with a function of a and /3. This is, on his turn, averaged over cr< giving the final result (Q"). 

A compact way for averaging over large times statistics is to introduce the Laplace transform of the PDF Vf = 
(exp[—s(5])>. It is indeed easy to recognize that its path integral coincides with a matrix element in the ’’Quantum 
mechanics” with Hamiltonian H = —d^ + s/3 exp {2r]) G [a exp ( 2 ? 7 ) /Pe^] /Pe^. The “space” variable rj = a{t') dt'. 

Using standard quantum mechanical notation, the expression for U/ can be presented as 


= (exp[-sQ])> = exp [-{T - to)/4] 


/ 


\ r 

(%) 

g-(T-to)Ao 

e-^/A = 




)/4$(r-to;??) 


r)—0 


(52) 


where $(t;? 7 ) satisfies (9* — H)^ = 0. The initial condition at t = 0 for the “wave-function” ^{t;ri) is exp (— 77 / 2 ). 
Noting that the potential part of the Hamiltonian vanishes ai rj ^ —00 and 4>(0; ry) does not, we obtain the asymptotic 
behavior of 4) at large times 


<^{t,7j = ln[?/Pe]) 


A/4 


^ 0 ( 4 /) 

VPey’ 


[d^ - sPG[ay‘^]] 4>o(j/) = 0. 


(53) 


The new variable y = (expry) /Pe has been introduced. The function $o(?y) should tend to unity for y ^ 0 and ^g/yty 
should vanish for y ^ 00 . It follows from (|^) that is simply the function 4)o calculated at ry = 1/Pe. 

The general way to attack ( |^ ) for an arbitrary form of pumping is to proceed as we have already done for ( |^ ) 
in Section Starting with a constant unit solution, the term with s in ( p^ is treated perturbatively. A series in s 
is then obtained U/ = 1-1- Cn{y)s'^ with the c„’s having an expression similar to (|H|) . The rz-th moment e" 

averaged over (T> coincides with c„(l/Pe), up to simple combinatorial factors. The final result is 

(e")> = (2rr-l)!!n! (54) 


where the coefficients a„ are given by 
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(55) 


2/2 


= / dy2G[y2] / dys 


OO 

J dyiG[y. 


y2n-2 


dy 2 n-l 


dy2nG[ylJ 


Equation (^, together with (H) allowing to calculate small times averages, gives the expression of all integer moments 
of the dissipation field for a general form of pumping. It follows from (54) that all these moments scale with the 
same power of Pe. This shows that the strong intermittency evidenced in the analysis of scalar differences in the 
convective range comes down into the dissipative range. The analysis of the n-dependence of the constants in (|5^ ) 
actually shows that the intermittency of the dissipation field is even stronger than for scalar differences. To this aim, 
it is convenient to restrict to a particular form of pumping, allowing to proceed with explicit calculations. Specifically, 
let us consider G*[x'^] = exp [—2a:]. The correlation function of the corresponding pumping has Fourier transform 
Xq = dy^/q'^exp (—l/g^). Equation ( ^ ) with G having the specific form G* is solved in terms of the Bessel function 
lo as 


If 


'P> = 

• c 


^ s/3/ae 
loWsPla) 


1 - 


Pe>i 


Pe lQ{^sf3/a) 


(56) 


The advantage with respect to the general case with arbitrary form of pumping is clearly that is now known 
explicitly and this will permit us to reconstruct the PDF of the dissipation field. 

Having averaged over cr>, we need now to take into account fluctuations at small times, i.e. average over (t<. 
The important remark is that, both in the general case (j^) and in (^^, we need to average quantities of the form 
^/^f(a/l3), with / arbitrary but having the property that it depends only on a/(3. For our purposes, it is then 
appropriate to introduce the random variable p, = a//3 and consider its distribution function weighted with 


V<[y\ 


f Va<exp(-S<) S f ({a<}) 

fVa< exp(-S<) 


(57) 


It is again convenient for averaging to introduce the Laplace transform of the PDF Vf = (exp (—s^))<. As in (|^), 
its expression can be presented as the following quantum mechanical matrix element: 


00 ! 

f e~^^'P^[y]dy = S{f]) 

exp 

—toH 

= 


d \ 

0 ' 



1 



(58) 


Here, satisfies (dt — 17)4' = 0, the Hamiltonian is H{ri;s) = —d^ + sexp(—2ry), the space variable y = J* a 

and the initial condition for the “wave function” is exp ( 77 / 2 ). The asymptotic behavior at large times t can 

be obtained as in (|^), noting that the potential part in H decreases at infinity and 'I'(0; y) does not. It follows 


4'(t;?7) 


e‘/‘^4'o(7?) = e^/'^exp [y/2 - Vsexp)-/;]] , 


(59) 


where 'I'o(? 7 ) satisfies (^H + 1/4^ 4/0 = 0 and behaves as exp ( 77 / 2 ) for large y's. Requiring exp (to/4) » I, we can 

plug the asymptotic expression (M) into ( ^^ and obtain Vf = exp [—y/sj. The PDF of y and the moments relevant 
for e" are easily restored as 


0'*’+ZCXD 


r<{y) = 


2Tri 




1 

27ri 




0+—i 


0+—2 


2^y?G 



j dyy-^+^/'^V<(y) = 2^^-^ 

0 


{n-iy. 



(60) 

(61) 


This expression can be used to calculate the moments appearing in (^). Note that the expression (|^), derived 
exploiting the time scales separation, coincides for the first moment (e) with the dominant order of the known 
solution derived in Section ^ : 
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e = 


1 f [x(0) - xjx)] dx 
7rPe J a;2 + Pe“^ 


(62) 


One step further can be made for the specific form of pumping G*[x^] = exp[—2a;], allowing to obtain the explicit 
solution (1^). From (|^), we obtain indeed the Laplace transform of the PDF for Q (averaged over both cr< and (t>) 
as 


OO OO 

Vs = {V>)<^l-^ [ ^ = 1 - [ ln[/o(2v^)]e"“dx. (63) 

Pe J V^PeJ 

0 XV// p 

The moments (e") can be immediately reconstructed, according to (^^, from the derivatives of ( |^ at s = 0. They 
read 


(e") 


(2n — l)!!n! 9" 
V^Pe ^ 


In 


1 

T{n+l/2)[T{n+l)f2^ 

1" \n[Jo{^/z)]dz 

JoiVz). 

* 27r2Pe j 

z=0 c 

ZU+I 


(64) 


where C is a close contour about z = 0 in the complex z- plane. The Cauchy integral representation ( j^ ) is useful for 
getting the large n asymptotic of (^) . The integral is indeed saddle around the square of the first zero of the Bessel 
function zq = Xq, Jo(xo) = 0. The saddle estimation for the integral is ~ Xg, thus giving 


ln[Pe(e")] 


n (31n[n] — 3 + ln[2xo]). 


(65) 


The whole expression (^3|) can actually be inverted. From (|4q), it follows indeed that the PDF of the 

dissipation e is given by 

0'*'+2OO 


0 0+—200 


= — / dzln 

Pe./ 




V2' 


1 


V? 


( 66 ) 


where qFp is the generalized hypergeometric function with the q parameters in the numerator and the p parameters 


in the denominator. The 5-function at the origin arising from the unit term in ( p^ has not been considered in 
The reason is that, as we shall see in a moment, the range of validity of ( |^ is e ^ Pe“^. 

It is indeed time to clarify the limits of validity of the calculations performed. To be concrete, we shall specif¬ 
ically consider the pumping G* leading to (^^. The remainder of the expansion over Pe“^ in (^^ is bounded by 

f3 \f^ /(2Pe^)- Contrary to (|^, we need now to consider averages oi p = a/(i with 

weight (3. The problem can again be reduced to the calculation of a quantum mechanical matrix element and it is 
found that the final result depends on tg as exp(2tg). The condition for the remainder to be subdominant with respect 
to the terms kept in (|^) is therefore that exp(2tg)/Pe should vanish as Pe ^ oo. On the other hand, the observation 
time must be much larger than unity in order to attain the stationary state, i.e. tg ^ 1. This condition is clearly 
compatible with the previous one, in the limit of large Peclet numbers, and gives the ordering 1 ^ to ^ ln[Pe]. 
The other delicate point is the criterium of applicability of (^ with respect to s. For the expansion in ( ^^ to be 
meaningful, the second term should be much smaller than unity. This shows that we should require s Pe"^. The 
expansion (^^ fails then to describe the large s tails of Vs and thus the smallest values of e in the respective PDF 
(the relation between large s and small e is direct since the decay of the generating function is relatively slow). 
Another simple approximation is however available for the high s limit. The trajectories contributing to (exp (—sQ)) 


at large s ^ Pe are clearly those where Q is small. From the definition (^) it follows that, for this to happen. 


W^{t) should remain small all the time. The quantity J W~‘^ is in this case 0(1) and the argument in G is 
small on account of the 1/Pe^ factor. For the trajectories relevant at high s, it is thus possible to approximate Q by 
G[0]/Pe^ f W^. The Laplace transform Vs reduces then to 
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Vs = ( exp 


[ dtW^it) 


Pe^ 


0 


- ^-T/4 


{ 


(5(7?) 


exp 


-TH{-r];sG[0]/Pe^) 


.-ri/2 


(67) 


where 7?(i) = f a{t) and H is the same as for (^). The matrix element (|^) actually coincides identically with 


when to by T, t] by —t], and s by sG[0]/Pe^ are replaced. Using (^, one can then easily get the final answer 


at s V> Pe 


exp 


\/sG[0] 


Pe 


( 68 ) 


Inverting (^), we can obtain the PDF of e at the small values e ^ Pe For larger values of e, the PDF follows from 
the general formula (|^. The following general behavior for the dissipation field PDF is thus obtained 


V2Pe 


V^^He) 


fo>'e«Pe2; 


1 - , 

111 v >„-2 


for Pe <e<l; 
, eo ^ 1, for e » 1, 


(69) 


exp[-(e/eo)^^^ 

Pi 


where algebraic prefactors have not been considered in region of exponential fall off. The PDF for scalar gradients 
uj = ^ye/K follows immediately from (|^ as 




1 


for oj 1/L; 

T’“(w) ^ for 1/L < w < I/Vk; 

^exp -(Iwly^K/eo)^/^ , for w » l/V^- 


(70) 


Moments of the dissipation (e“), with a > 1/2, are all proportional to 1/Pe, in agreement with (54) valid for arbitrary 
pumping. Moreover, we learn from (|^ ) that the linear scaling found for scalar differences low-order moments comes 
down to the dissipative range. Moments (e“) with —l/2<a<l/2 scale indeed as Pe“^“. Moments with a < —1/2 do 
not exist. ?,From (70), it also follows that the same tendency observed for scalar differences PDF to develop a cusped 
structure at the origin is present. An important difference arises for the tails. The comparison between and 
( ^ ) for gradients and scalar differences indicate indeed that the former decrease much more slowly. Such behavior is 
physically understood in terms of fluctuations of the dissipative scale as follows. Gradients can be thought of as scalar 
differences evaluated at the dissipative scale. This scale is however a dynamical quantity and fluctuates (strongly in 
a problem with inverse transfer as the one here). The statistics of these fluctuations are actually related precisely 
to those of the variable /r, introduced in (p^. The fluctuations of e" are therefore the product of those of scalar 
differences and those of the dissipative scale. The latter do not depend on the dimensional parameters of the problem, 
i.e. they are universal as follows from (|^), but grow factorially with the order n. This factorial is precisely what 
shifts the dominant term in (^5|) from 2nhin to 3nlnn. It follows then that the slower decay of the tails for gradients 
with respect to those for scalar differences is indeed due to dynamical fluctuations of the dissipative scale. 


7. CONCLUSIONS AND DISCUSSION 

We analyzed in the present paper statistics of the passive scalar advected by a one-dimensional, smooth and fast- 
correlated in time velocity field. In spite of its extreme simplicity, the dynamics presents very interesting and surprising 
behaviors. We expect that many of them are quite generic and it may be interesting to look for more realistic scalar 
(and generally turbulent) problems where the new effects discussed here could be of relevance. Thus, hereafter we 
recall and briefly comment the major properties of the model discussed in the paper and address, in parallel, the 
respective questions for future studies. 

The first major property of the dynamics is the inverse cascade of the scalar. The inverse cascade is a consequence 
of compressibility, but in a subtle way: in we consider indeed a generalized smooth d—dimensional model having 
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the degree of compressibility as a free parameter. A transition between inverse and direct cascades is observed there : 
if d > 4, the cascade is always direct, independently of the degree of compressibility; if the latter is small enough, the 
cascade is direct again; otherwise, it is inverse. It might also be interesting to look from the same dynamical point 
of view used here at the critical dimension appearing in the passive vector problem discussed in |25 and at the role 
of compressibility for the direction of energy transfer in other turbulent situations, e.g. 2D Navier-Stokes and MHD 
turbulence. 

Dynamically, the inverse cascade is associated with the fact that the Lyapunov exponent for Lagrangian trajectories 
is negative. This implies that, in the case of inverse cascade, contraction of Lagrangian trajectories is typical and 
stretchings are relatively rare events. This is precisely the opposite of the canonical picture associated with the direct 
cascade. Rare trajectories (contracting or stretching if the Lyapunov exponent is positive or negative, respectively) are 
responsible for the intermittent part of the problem. This definetely confirms the picture that contractive trajectories 
play a crucial role for intermittency in incompressible turbulent transport. 

Negative Lyapunov exponent leads to the third observation emerging from the analytical study of the model: 
Gaussianity of the scalar establishes at scales larger than the scale of the pumping, while strong intermittency is 
present at small scales. Small-scales intermittency found here is of the “Burgers” kind, i.e. all scalar differences 
integer moments scale linearly ^ r/L. This result gives a definite ground to speculate that, generally, the stronger is 
the degree of compressibility (the role of contracting trajectories is then enhanced), the more intermittency downscales 
of the pumping is observed. Note in the respect of Gaussianity that the inverse cascade of energy in 2D Navier-Stokes 
turbulence is also Gaussian according to the numerics |3|]. It rises yet another conjecture of a certain generality of 
the Gaussianity feature for a great variety of direct cascades. 

The equation for scalar differences PDF, derived directly from the dynamics in Section shows that no dissipative 
anomaly is present in the model. In the generalized one-dimensional model introduced in |^, 7 = I is a natural 
threshold for the inverse cascade. At 7 > I the cascade is direct, while it is inverse at 7 < 1. Most probably, it 
means that there is no dissipative anomaly (equation for the PDF of the scalar difference could be closed) for 7 < 1. 
However, it is not yet clear how to derive this consistently (as it is derived here for the smooth limit 7 = 0 ). Note 
that the general question of the role of the dissipative anomaly in the passive scalar physics is not yet resolved in any 
sense (some suggestions on this point may be found in [|^,|^,^). One more unresolved question is how does the 
dissipative anomaly (if any) affects the anomalous scaling behavior of structure functions of integer (and generally all 
the) orders. 

An important result for the issue of convection-diffusion interplay, is the fact that the effective dissipative scale, 
which sets the crossover scale between scalar differences 2 n-th moment in the convective range and dissipation field n- 
th moment, is a strongly fluctuating quantity, growing factorially with n. To make this statement, we have analytically 
calculated in Section ^ the dissipation field PDF exploiting a new scale-separation procedure, which is also used in 
Appendix B to describe the long-time dynamics of the pair correlation function of the scalar gradients. The procedure 
is likely to be a general tool for resolving the problem in other situations. Particularly, a slight (matrix) modification 
of the method should be relevant for the Kraichnan d—dimensional incompressible smooth model. 

Finally, the eddy-diffusivity operator is Hermitian for Kraichnan model, where random velocity is incompressible 
and fast. The Hermiticity is lost, when compressible flow is considered. Generally, it might be interesting to look 
at scalar transport (and particularly at flow with finite correlation times) from the point of view of non-Hermiticity. 
This issue is the subject of a very recent interest in the field of disordered systems The specific object 

of interest for scalar transport is the resolvent TZ(t]x,y), describing the probability to find a Lagrangian separation 
equal to x at time t, having initially been equal to y. One may then want to study the distribution in the complex 
plane of the poles of its Laplace transform. A width of the poles domain along the imaginary direction might be a 
characteristic of the trapping degree. 
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APPENDIX A: DYNAMICAL FORMULATION BY FIELD FORMALISM 


Let us consider the gradient field uj{t,x) satisfying the equation of motion with the velocity u{t,x) having 
Gaussian statistics and correlation function < u{t,x) u{t',x') >= V{t — t',x — x'). The generating functional T (A) 
for its simultaneous correlation functions can be written in the form of functional integral (see 


iF (A) = J exp i (p|x"|p) — i (u|U J dtdxp (^dtuj + dxiuLo) — nd'^co) + J dxX{x) lo {T, x) 

(Al) 

Here, [u\V~^\u) denotes the diagonal matrix element of the inverse operator to V{t — t',x — x') in the Hilbert space of 
the functions u (t, x) and the pumping y" is the second spatial derivative of the forcing correlation function y, defined 
in The retarded regularization of the time derivative, (dtuj)^ = {tUn — ijJn-i), does not produce any nontrivial 
Jacobian after derivation of from the equation of motion. A 0 is the temporal slicing and (cc) = to (nA, x). 
Performing the Gaussian integration over the field u (t, x) we arrive at 

.F (A) = J UpUaj exp | —+ 2 (fIx^Ip) + * J dtdxp (dtUJ — kOIuj) + J dccA (a;) w (T, a;) | . (A 2 ) 

For the case of J-correlated in time velocity fields considered in this paper, we replace the general V{t — t',x — x') by 
specific 6 {t — t')V{x — x'). The explicit version of & reads then as 


5 = 


T (A) = J VpVuj exp + J (^) ^ 

T 

- y dt J dxdx'{uidxp) {t,x)V {x — x') (uidx'p) {t,x') — J dxdx'p {t,x) x'‘ 
0 

T 

—i J dtdxp [dtijJ — nd^oj) ■ 


(A3) 


X — X 


p{t,x') 


(A4) 


Despite of the explicit presence of the U-term in (A4), the dynamical part of the action S is covariant with respect 
to Galilean transformations. This leads to the important Ward identity: 


J 'Dp'Du! exp + J (^) ^ n(/ dxuj{tj,x)dxp{tj,x)^ = 0 , 

where set {tj < T,j = 1, 2,...} is arbitrary. To prove it, let us consider the change of variables 

( } \ ( } ' 

p {t, x) — > p \ X + v{T)dT, t \ , uj {t,x) — > a; a; + / v{T)dT, t 


(A5) 


(A 6 ) 


in the functional integral (A3) with v (t) being unconstrained. The transformation does not change the source term 
and it has unit Jacobian. The variation 5S of the action S under the transformation reads as 


T 

SS =—i J dtv{t) dxuj {t,x) dxp{t,x)^ . 


(A7) 


The change of variables does not change the value of the integral. Thus, all the functional derivatives of (A2) over 
v{t) are equal to zero. Taking into account (A7), one finally arrives at the Ward identity (A5). 

All written above is applicable for any short correlated velocity field. Now let us narrow consideration to the case 
of smooth velocity field (® with 7 = 0 . One gets 
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JF (A) = J 'Dp'Du! exp ^—S + J dxX{x) to {T,x)^ , 

T 

— J dxdx'(ujdxp) {t, x) {x — x')^ (ujdx'p) (t, x') — J dxdx'p {t, x) x" 

0 

T 

—i J dtdxp {ptLo — , 


S=- dt 


X — X 


p{t,x') 


(A8) 

(A9) 

(AlO) 


where the Vq term was dropped due to (A5). The goal is to reduce the problem to averaging of a functional of 
one random in time scalar field, analogously to the random matrix description used in two- dimensional and 

generally d-dimensional case. However, in the present case, compressibility calls for more evaluations. Let us perform 
the following change of variables in (A8): 

{p{t,x), uj{t,x)} i^p ^t,x + i J dr J dy y'^dyp{t,y)uj{t,y)^ , uj x + i J dr J dy y'^dyp{t,y)u}{t,y) 

(All) 

which can be considered as the field-dependent and time-dependent homogeneous spatial translation leaving the source 
term intact. Analytic continuation of the fields which makes the spatial arguments to be real is assumed. Action gets 
the following from being rewritten in the new variables 


S = \jdt 


2 1 ^ 

dxx{ujdxp) {t,x)^ — J dxdx'p{t,x)x" J dtdxp {^dtoj — kO^uj^ . (A 12 ) 


The transformation (All) has non-trivial Jacobian 

'Dp'Doj —> 'Dp'Duj J (p, uS ), 


(A13) 


which depends on regularization. The regularization is fixed by the requirement for the temporal <5— function from the 
correlation function of velocities to appear as a result of narrowing of an even function of temporal argument. Besides, 
the correctly regularized action (A12) should reproduce the respective correlation function in the limit T —> 0: 


N 


Pn {x) — >pnix + - A / dy y'^dyPn {y) ujn (y) -I- iA ^ dy y'^dyPm (y) uJm (y) , 


respectively for Un (x). The Jacobian ^7 (p, a;)can be computed easily 


(p, w) = exp z J dtdx xudxP^ ■ 


(A14) 


(A15) 


T ^ 

The exponential of the term f dt (f dxxLudxp) in the action (Al2|) multiplied by the Jacobian (A15) can be rep- 

0 


resented by means of averaging of the exponential of i f dta f dxxwdxP over the auxiliary field The Gaussian 

0 

measure of averaging over a is defined in (pi[). It results in the following representation for the generating functional 


-5.--(/|x-V) 


T {\) = J 'Dp'Dui'Df'Daexp 

T 

i J dtdxp (ptOJ — adx{xuj) - UJ — dxf) + J dxX{x)uj {T,x) 


xexp 


(A16) 

(A17) 
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Integration over p (t, x) gives the reduced equation of motion: 


dtoj - ada:{xuj) - ndloj = dxf, (A18) 

averaged with respect to / (t, x) and a (t) with the weights given by the first line of ( |A17| ). It should be stressed that 
this equivalence holds only at the level of simultaneous correlation functions. 


APPENDIX B: LONG-TIME BEHAVIOR OF THE PAIR CORRELATION FUNCTION OF GRADIENTS. 


In the present Appendix, we will find the long-time asymptotic of the pair correlation function of the scalar gradients 
n (t, x) which obeys (i- We will model here the infrared cut-off of the velocity correlations by 


The Laplace transform of 


S (x) = x^, X < Lu, S (x) = LI, X > Lu- 
has the form 

sfts — dl {S {x) + 2 k) ils + x" {^) /s = 0. 


(Bl) 


(B2) 


The solution of (|B2| ) can be written in terms of the Green function Qs{x^x') (which is taken to be even in x, 
Gs {x, -x') = Gs {x, x')) as 


OO 

fls {x) = - - J Gs {x, x') x" (x'^ dx , sGs - dl {S (a;) + 2k) Gs = S (^x - x ^ . 


(B3) 


The desired long-time asymptotic corresponds to the smallest s and our aim is to find the singular structure of Gls (x) 
at small s and x ^ ^/K . In what follows, we will then neglect all the logarithmic contributions like s In and s In a: 
in comparison with 1. We will also neglect the smallest x' {x' <C Vk) contribution into the integral going to 
zero in the limit of small diffusivity. Therefore Gs (x,x') should be studied at ^/k x' and in all the allowed 

domains with respect to x. 

Making the substitution Gs {x, x') = g {x, x') /[S' (a:) -I- 2k], one gets 


S (a:) -I- 2k 


9 -dlg = 5{x- x^ . 


Solution of (^) in the domain of the largest x > Lu > x' ^ can thus be written as 


■\/s 

g {x,x') = Aexp ( (a; - L„) 

Lu 


(B4) 


(B5) 


Note however that only the first two terms of the ( p35| ) expansion in s will be required for matching hereafter. If x 
is smaller than Lu, but is still larger than a separation scale xq, one can omit the hrst term on the rhs of (B4) and 
write the following general solution of (B4) in the intermediate domain, at Lu > a; 3> ^/k, 


g{x,x') = 


Fi + xF 2 , X < x', 
Di + xD2, X > x' . 


(B6) 


Here, the constants Fi, F 2 , Di, D 2 should be fixed via inner matching at a: = a;' and outer ones at Lu and xq. The 
value of Xq, restricted by ^/k xq ^ \fi^ls, defines the upper bound for the scales at which one cannot omit the first 
term of (B4) anymore. However, the point is that at x < xo < x' an actual dependance of g on xo has disappeared 
and we should not locate xo explicitly to get the asymptotic behavior of g{x,x'). Indeed, integrating we get 


X 

9x9 {x) = s J 


9 jz) dz 
S(x) +2 k ’ 


(B7) 


where the argument x' is omitted and it is already taken into account that g{x) is even: 8x9 (x) U=o = 0- (B7) shows 
that the x—derivative changes sufficiently at the diffusive scale, being nonetheless small by the absolute value there. 
It means that the function g (x) at x xq is a constant in the leading (zero) order in s. The subleading (first order 
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in s) term is then fixed by (B7) taken at x = Xq. Alltogether it gives the following expression for g at x < x' with the 
desirable (for the forthcoming matching) accuracy 


g~ go 1 + sa; 


dz 


S{x) + 2k 


50 1 + s 


2V^) 


where, evaluating the integral in (^^ ), we used that ^/k <C xq. Matching (Bf.and (^), we get 

at X = x' Di+ x'D2 = Fi+ x'F2, F2 - D2 = 1 , 

at X = Fu Di + L^D2 — A, D2 = —A, 

at X = Xq + X0F2 = go, F2 = sgo^. 


Keeping s to be very small, while the ratios sjy/K and y/sjLu are finite, one derives from (B9) 


Fi = go = 


1 + x' ^/sj Lu 
s-k/ {2rd) + ^/s/L, 


F2 = 


SgpTT 

2V^’ 


^ . 2y^ + x'sTr y/s . 

Di=A= --=^—, D2 = -^A. 

STT + 2',/2ksILu 


It is worth noting that the explicit value of x^ does not enter (BIO), and the scale separation procedure 
justihed. Substituting (BIO) into (^) (x' < L^), one gets 


(x) — 


S (x) + 2k 

and coming to the inverse Laplace-transform 

1 

Vt (t, x) = 


X{x) x(0) 


S + ay/s 


a = 


2-\/2k 

'kL,, 


S (x) + 2k 


0 


p^ + l 


The expression (|B12| ) states that, at the times less than 7^^ ^ the correlation function coincides 
quasi-stationary limit given by ( 0 )- At the largest times, (t, x) approaches the asymptotic value 


1 

S (x) + 2k 


X{x) 


diffusively, x(x) — (x) + 2k\ = C ^ x (0) (a^Trt) 
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